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ABSTRACT 

Within a cosmological hydrodynamical simulation, we form a disc galaxy with sub- 
components which can be assigned to a thin stellar disc, thick disk, and a low mass 
stellar halo via a chemical decomposition. The thin and thick disc populations so se- 
lected are distinct in their ages, kinematics, and metallicities. Thin disc stars are young 
(<6.6 Gyr), possess low velocity dispersion (ctu,v,w = 41, 31, 25 kms -1 ), high [Fe/H], 
and low [O/Fe]. Conversely, the thick disc stars are old (6.6<age<9.8 Gyrs), lag the 
thin disc by ~21 km/s, possess higher velocity dispersion (ou,v,w = ^9, 44, 35 kms -1 ), 
and have relatively low [Fe/H] and high [O/Fe]. The halo component comprises less 
than 4% of stars in the "solar annulus" of the simulation, has low metallicity, a veloc- 
ity ellipsoid defined by (cttj,v,w — 62,46,45 kms -1 ) and is formed primarily in-situ 
during an early merger epoch. Gas-rich mergers during this epoch play a major role 
in fuelling the formation of the old disc stars (the thick disc) . We demonstrate that 
this is consistent with studies which show that cold accretion is the main source of 
a disc galaxy's baryons. Our simulation initially forms a relatively short (scalelength 
~T.7kpc at z=l) and kinematically hot disc, primarily from gas accreted during the 
galaxy's merger epoch. Far from being a competing formation scenario, we show that 
migration is crucial for reconciling the short, hot, discs which form at high rcdshift in 
ACDM, with the properties of the thick disc at z—0. The thick disc, as defined by its 
abundances maintains its relatively short scale-length at z = (2.31 kpc) compared 
with the total disc scale-length of 2.73 kpc. The inside-out nature of disc growth is im- 
printed the evolution of abundances such that the metal poor a-young population has 
a larger scale- length (4.07 kpc) than the more chemically evolved metal rich a-young 
population (2.74 kpc). 

Key words: galaxies: evolution - galaxies: abundances - methods: numerical 



1 INTRODUCTION 

The fitting of the vertical distribution by two exponentials 
has motivated separation of the Milky Way's disc popula- 
tion into 'thick' and 'thin' stellar disc components. Such 
classification has dominated the liter ature of Galactic struc- 
ture for the past three decad es (e.g. iGilmore fc Reidlll983l; 
Juric et al]|2008l . although see llvezic et al.ll2008l : lBovv et all 
2012j h The oldest disc stars in the Milky Way are kinemat- 
ically hot, show an asymmetric drift (i.e. rotational lag), 
and have lower metallicities ([Fe/H]) and higher a-to-iron 



abundance ratios, with respect to the younger disc stars. 
Understanding the origin of the stars distributed at the 
highest vertical scaleheights within the disc - the classi- 
cal 'thick disc' - is an important aspect in our attempts 
to understand galaxy formation and evolution throughout 
the Universe, as such thick discs appear ubiquitous (e.g . 
iDalcanton fc Bernstein] 12002; ; lYoachim fc Dalcantor] 120051 ). 
and contain the oldest disc stars, born during the earliest 
formation of discs. 

Several models for the formation of such thick stellar 
discs have been proposed, each of which are able to repro- 
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duce one or more of their empirical properties. Despite these 
partial successes, a singular, universally accepted model has 
not yet emerged. Such a model needs to answer two funda- 
mental points: (i) the origin of the oldest disc stars, and (ii) 
how these stars attained their present day properties. 

The heating of a pre-existing thin disc due to 
the accretion of satellite galaxies (jQuinn et all ll993T ) 
appears consistent with the hierarchical structure for- 
mation of ACDM. The signatures of such events are 
well studied and indeed are able to match verti- 
cal profiles of the thick disc (Havashi & Chiba 2006: 
Kazantzidis et al.ll2008l: IVillalobos fc Helmill200Sl : iQu et all 
201 ll : iBekki fc Tsuiimotol |201lT) including several fine- 
detailed features I Qu et al.l l201ll ). Such models assume a 



pre-existing and relatively thin disc. Another model pro- 
poses that stars accreted by satellites can be dragged 
into the plane of a p re-existing disc by dynamical fric- 
tion l|Abadi et al.ll2003h . Recently, it has been shown that 
the radial migration of old, a-enhanced stars from the in- 
ner region of the disc to the solar nei ghbourhood can ex- 
plain many properties of th e thick disc ( Schonrich fc Binnevl 



l2009l: iLoebman et alj|201ll). T he "po pping" of star clusters 
[ Kroupall2002c 1 Assmann et al.lfeoill ) or the effects of mas- 
sive clump formation in unstable gas-rich discs in the early 



universe (TNoguchill 19991; iBournaud et alj|2007l ; lAgertz et alj 
120091 : ICeverino et al.l 120101 ) may also result in thickening 
of discs; such scenarios link thickening more to formation 
processes, rather than to disc heating, per se. Stars may 
also have been born relatively thick as gas settles into a 
disc configuration du ring an early period of gas-rich merg- 
ers l|Brook et al.ll2004bl ). 

It is important to stress that these models are not 
necessarily mutually exclusive. Indeed, some do not try to 
explain the formation of the oldest stars, but rather how 
they were heated, or how they migrated to the solar re- 
gion. One way to comprehensively probe thick disc forma- 
tion is using high-resolution, hydrodynamical, cosmological 
galaxy forma tion simulations of Milky Way mass galaxies 
( Ab adi et al.l 120031: IStinson et al. 2010; Scannapicco et al. 
20101: iKobavashi fc Nakasatd l201ll: iTissera et al.l l2oilT 
Martig et al.ll2012l : iDome ncch-Mo ral et al.ll2012T ). These at- 



tempts, though, have been hampered by their inability to 
match several important properties of the Milky Way. In 
this context, the main problems are (i) the continued for- 
mation of overly massive and metal rich stellar halos, which 
makes separation of thick disc and halo stars more diffi- 
cult than in the Milky Way where metallicity can be used 
as an added discriminant, and (ii) the lack of resolution. 
The first problem has been alleviated by improved feedback 



llBrook et al.ll2004al: lOkamoto et~al] 120051 : iGovernato et all 
120071 : iGuedes et al.[|201ll) bu t con t inues to plague even re- 



cent simulations jSales et al.ll201ll : IScannapieco et al1l201ll : 
iMartig et al.ll2012l ). The latter problem is more subtle than 
just the size of gravitation softening, but relates to the rel- 
atively low density and high temperature at which star for- 
mation occurs in low-res olution simulations , which prevents 
thin discs from forming l|House et al.ll20l"lh . 

In this study, we take a complementary approach by 
examining a disc galaxy simulation of lower mass than the 
Milky Way, providing us with improved spatial resolution 
and allowing stars to form at higher densities and lower 
temperatures. We analyse the simulation first presented in 



iBrook et al.l (|2012al . B12, hereafter) late-type barred galaxy 
with M„ ir =1.9xl0 11 M , realised as part of the MaGICC 
(Making Galaxies In a Cosmological Context) program. To 
re-iterate, this simulation is of a lower mass than that of 
the Milky Way, and therefore our analysis must be read 
within this context. That said, we note that the simulated 
galaxy examined here shares other properties with the Milky 
Way that we feel make it relevant to analyse in relation 
to putative theories of thick disc formation. Although the 
Milky Way is the only galaxy for which we have detailed 
abundances and kinematic information, vertically-extended 
'thick' discs appear ubiquitous in nature. Any comparison 
we make with the Milky Way is not expected to match ex- 
actly, but is predicated upon the fact that the broad trends 
of the Milky Way chemo-dynamics should be generic for 
relatively massive disc galaxies which have quiescent merger 
histories, in the sense that the last major merger was early in 
the galaxy's evolution. There is evidence that the Milky Way 
had such a quiescent merging history (e.g. lHammer et al.l 
120071) . However, if these assumptions were to be shown to 
be false, our model would clearly unravel. 

The second feature of our study is our chemical de- 
composition of stars into different galactic sub-components 
based on abundance patterns. The ability to separate com- 
ponents based purely on the ir chemical "tags" is at the 
heart of Galactic Archeaology (|Freeman fc Bland-Hawthornl 
2002). A multitude of studies has shown that thin and 
thick disc stars with similar iron abundances ([Fe/H]) 



differ in their abundances of a elements (e.g. iMai ewski 



1 1993 1; | Fuhrmannlll99Sl; IChiba fc Beers! [20001: iBensbv et al 
hood ; ISoubiran et alj|2003l ; IWvse et al.l 120061 : iReddv et al 



2006; iFuhrmannl 120081 : iRuchti et al.ll201Ct). Using a com- 



pilation of observational data, Navarro et al.l (|201ll ) made 
simple cuts in the [a/Fe]-[Fe/H] plane, in order to sepa- 
rate the Milky Way's thick and thin discs. We adopt a 
similar decomposit i on in t he simulated [0/Fe]-[Fe/H] plane. 
iBensbv fc Feltzind l|201ll ) argue that the large overlap in 
kinematics between thick and thin disc populations makes 
classification based on age better than classification based 
on kinematics. The relative difficulty of ascertaining stel- 
lar ages means that abundances, in conj unction with abun- 
dance ratios, may be the best "proxy" . iBovv et all (l201ll ) 
make detailed cuts in abundance space to probe thick disc 
properties. 

Thirdly, our study proposes more detailed categorisa- 
tions of the origin of stars that comprise the different compo- 
nents, going beyond accreted vs in situ. In particular, we con- 
cern ourselves with whether gas which subsequently forms 
stars has been involved in a "gas-rich merger", as well as 
whether it was smoothly accreted or accreted directly from 
proto-galaxies/satellites. This provides a more complete per- 
spective of the relative importance of gas-rich mergers com- 
pared with smooth accretion. 

In Section [2 we review our code and initial conditions, 
and outline the basic pr operties of the r esulti ng simulation. 
In Section \2. 21 we follow H avarro et al.1 l|201ll) and separate 
our (analogous) solar neighbourhood stars on the basis of 
their abundances, and find that such classification results 
in components that could reasonably be referred to as a 
thin disc, thick disc, and halo/metal weak thick disc, shar- 
ing some gross characteristics to those in the Milky Way. 
In Section [4] we trace the origin of the stars in these pop- 
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Figure 1. The evolution of the simulated galaxy shown in B- 
band surface brightness maps, as viewed both face-on (left panels) 
and edge-on (right panels). The redshifts range from z=2.7 to 
2=0, and are noted in the bottom right corner of the edge-on 
panels. Each panel is 40x40 kpc (physical units). The effects of 
dust reprocessing, as implemented with Sunrise, are included. 



ulations, both in terms of where they formed (in situ or 
accreted, and where in the disc that they formed), their for- 
mation dispersion, and the origin of the gas from which they 
formed. We find that thick and thin disc stars form in situ 
and that cold accretion and gas-rich mergers supply thick 
disc gas. We find in Section [5] that migration also plays a 
role in forming the simulated thick disc. In Section [7J we 
discuss the consistency of these results and the fact that 
most disc stars form from smooth accretion and highlight 
the role of migration in shaping our thick disc population. 



2 THE SIMULATION 

2.1 Code and Initial Conditions 

We analyse the simulation described in B12. To quickly 
review, this is a cosmological 'zoom' simulation, derived 
from the initial conditions associated with galaxy gl5784 
from the McMaster Unbiased Galaxy Simulations (MUGS, 
IStinson et afll20ld ). using a 25 h 1 Mpc parent cube. The 
virial mass of the simulated galaxy is M V i r — 1.94 x 10 Mq . 
Further detailed chemical modelling (see below) has been 
incorporated into the code, and the simulation re-run: the 
stochastic nature of our star formation recipes mean that 
small differences exist between the current version of the 
simulated galaxy and the one analysed in B12, but essen- 
tial properties such as rotation curve shape, star formation 
history, magnitude, colour, and disc scalelength are not al- 
tered significantly. In B12, we showed that secular processes, 
largely driven by a strong bar, resulted in the formation of 
a bulge with B/T=0.21 at z=0. In the current simulation, 
a strong bar does not form and hence, neither does such a 
significant bulge. The current simulation has B/T=0.12 at 
z=0. 

The simulation is evolved with t he smoothed particl e 
hydrodynamics (SPH) code gasoline |Wadslev et al.l2004T l. 
GASOLINE employs coolin g based on the r adiative transfer 
treatment within CLOUDY l|Shen et al.ll2010t ). including cool- 
ing due to H, He, a variety of metal lines, in the presence 
of an external UV radiation field. A pressure floor is added 
^Robertson fc Kravtsovl2008l ) , and a maximum density limit 
is imposed by setting a minimum SPH smoothing length, 
ham to i of the gravitational softening to ensure that gas 
resolves the Jeans mass and does not artificially fragment. 
We have ~ 4 million resolution elements (gas+dark mat- 
ter+star) within the virial radius at z — 0, with mean stel- 
lar particle mass of 7300 Mq and a gravitational smoothing 
length of 155 pc. 

When gas reaches cool (T< 10,000 K) temperatures in 
a dense (nth>9.3 cm -3 ) environment, it becomes eligible to 
form stars. This value for nth is the maximum density gas 
can reach using gravity. Gas is converted to stars according 
to the equation 



AAf* 



tdy 



(1) 



where AM* is the mass of the star particle formed, At the 
timestep between star formation events, M gas the mass of 
the gas particle and tdyn the gas particle's dynamical time, 
c* is the efficiency of star formation; in other words, the 
fraction of gas that will be converted into stars during tdyn ■ 
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Figure 2. [0/Fe]-[Fe-H] of simulation stars in the region 7 < Rxy < 8kpc, \z\ < 0.5kpc. Two "evolutionary tracks" are evident, 
loosely delineated by the dashed line: [O/Fe]=-0.24x [Fc/H]— 0.16. This simple delineation defines what we label 'thick' and 'thin' disc 
populations in this work. We then overplot the means of the thick (red diamonds) and thin (blue triangles) discs [O/Fe] as a function of 
[Fe/H], for the populations as defined by their ages (see text). Above (on the right) is plotted (on a log scale) the distribution of [Fe/H] 
([O/Fe]) in the 'solar annulus' for thin disc (blue line), and thick disc (red line) stars, based on abundances (solid) and age (dashed). 



Two types of energetic feedback are considered: su- 
pernovae and stellar radiation. Supernova feedback is im- 
plemented us i ng th e blast-wave formalism described in 
IStinson et al.l |200q ) and each supernova is assumed to de- 
posit 10 51 erg of energy into the surrounding interstellar 
medium (ISM). Since this gas is dense, the energy would 
be quickly radiated away due to the efficient cooling in 
SPH. For this reason, cooling is disabled for particles in- 
side the blast region. We model the luminosity of stars us- 
ing the iTorresI (|2010h mass-luminosity relationship. These 
photons do not cou ple efficiently with the surrounding ISM 
|Frever et alj|2003 l. We thus want to couple only a small 
fraction of this energy to the surrounding gas in the simu- 
lation. To mimic this highly inefficient energy coupling, we 
firstly inject 10% of the energy as thermal energy in the sur- 
rounding gas, and cooling is not turned off for this form of 
energy input. It is well established that such thermal energy 
injection is highly inefficient at the spatial and temporal 
resolution o f the type of cosmological simulations used here 
jKatdl 19921 ). This is primarily due to the characteristic cool- 
ing timescales in the star forming regions being lower than 
the timestep of the simulations. 

Metals are ejected from type II supernovae (SNell), 
type la supernovae (SNela), and the stellar winds driven 
from asymptotic giant branch (AGB) stars. The ejected 
mass and metals are distributed to the nearest neighbour ga s 
particles using the smoothing kernel (|Stinson et al.l 12006). 
and w e employed standard yie lds from the literature for 
SNel l (jWooslev fc Weaverll 19951 ) and SNela jNomoto et all 
1997). Metal diffusion is also included, such that unresolved 
turbulent mixing is treated as a shear-dependent diffusion 



term (Shen et al. 20l"oT). This allows proximate gas particles 
to mix their metals. Metal cooling is calculated based on the 
diffused metals. 



2.2 Evolution and Properties 

We emphasise again that this is not a simulated Milky 
Way galaxy. Our simulated galaxy is a less massive late- 
type galaxy, with rotation velocity of ~140 km/s, and a 
total stellar mass within the virial radius of ~8xl0 9 Mq. 
The galaxy has an absolute magnitude of Mb, v,i=— 19.45, 

— 19.94, —20.71, and sits in the blue sequence of the colour- 
magnitude diagram of observed galaxies. 

We plot the morphological evolution of the central 
galaxy in the B-band in Figure [T] face-on (left panels) and 
edge-on (right panels) from 2=2.7 to 2=0 (bottom to top 
in the figure). Each panel is 40x40 kpc (physical units). 
The major merging activity of the galaxy occurs between 
2=2.7 and 2=1.7, and we refer to this period as the "merger 
epoch". Three significant mergers occur during this period 

- firstly, a 3:1 major merger, followed by mergers with mass 
ratios of 10:1 and 20:1, with the last vestiges of a merging 
satellite still apparent at 2=1.5. In B12, we showed that over 
90% of the gas which cools to the bulge region during the 
merger epoch is blown out to the hot halo or beyond, rather 
than forming bulge stars. The galaxy does not form a signifi- 
cant bulge during the merger epoch and has an almost pure 
exponential profile at 2=1 (Sersic index n=l.l), at which 
time it has a scalelength of 1.7 kpc. The scale-length of the 
disc component at 2=0 is 2.7 kpc, as measured in the /-band. 



Disc Galaxy Decomposition 5 



Redshift 
0.5 1 1.5 2 3 
0.03 P ' 1 1 1 r 




2 4 6 8 10 12 
Stellar Age (Gyrs) 

Figure 3. The star formation history (SFH) of stars in the solar 
annulus, defined as the region 7 < R xy < 8 and \z\ < 0.5 kpc 
(black line). The SFH of the thin disc (blue line), thick disc (red 
line), and halo (yellow line), as defined by abundances in the 
[0/Fe]-[Fe-H] plane, are overplotted. 

3 RESULTS 
3.1 Components 

In Fig we plot [O/Fe] vs [Fe/H] in the 'solar annulus', 
taken as the region 7 < R xy < 8 kpc and \z\ < 0.5kpc, 
where the disc is in the xy plane Q We analyse here oxygen 
as a representative a element which is synthesised predomi- 
nantly in Type II SNe, and which is an element whose syn- 
thesis is relatively well studied. An interesting feature of this 
plot is that there appears to be two trajectories, one stretch- 
ing from ([Fe/H], [O/Fe] )~ (-1.1,0.11) to ( [Fe/H], [O/Fe] )~ 
(-0.82,0.05), and the other running from ( [Fe/H], [O/Fe] )~ 
(-0.85,0.035) to ( [Fe/H], [O/Fe] )~ (-0.5,-0.02), offset by 
~0. 01-0. 05 dex in [O/Fe]. We explore these trajectories by 
separating them into two populations simply using a straight 
line with equation [O/Fe]=-0.24x [Fe/H] -0.16 (the dashed 
line in Fig [2j . 

The ~0. 01-0. 05 dex vertical offset between the two 
'tracks' appears to be the result of the reduction in the star 
formation rate (SFR) ~6.6 Gyrs ago, when [Fe/H]~— 0.85 
(see Figs. [7]&[3} . In the case of the Milky Way, the tracks are 
offset vertically by ~0.2 dex, rather than the ~0. 01-0. 05 dex 
seen here. Yet the existence of these two tracks and the na- 
ture of their offset motivates us to associate the populations 
with the thick (above the line) and thin (below the line) 
discs. 

The larger offset in the Milky Way (0.2 dex) has been 
interpreted as indicating a hiatus in star formation in the 
Milky Way, along with an initiation of metal-poor gas in- 
flow associated with the formation of the thin disc (e.g. 
IChiappini et al.ll200ll ; iFenner fc Gibsonll2003h . Increased in- 
flows of metal-poor gas are not seen ^6.6 Gyrs ago in this 
simulation, and we see a star formation reduction rather 
than hiatus. This could explain the difference with the 
dual infall models in the size of the offset. Further simu- 
lations with more massive galaxies will be needed to explore 

1 We did similar analyses for various regions between 5 < R xy < 
9 kpc and \z\ < 2 kpc; none of the trends discussed are dependent 
upon the precise region chosen as representative of the disc. 
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Figure 4. The distribution of rotational velocities in the solar 
annulus for the thin disc (blue line) and thick disc (red line) 
stars, based on abundances (solid) and ages (dashed). Note the 
log scale. 



whether a star formation reduction coupled with the gas 
cycle produced in our simulation can result in larger off- 
sets in the tracks, or whether further processes can be in- 
voked tha t would lead to a h i atus in star formation. We 
note that ISchonrich fc Binnevl (2009) attain bimodal dis- 
tributions in [O/Fe] as a consequence of their assumptions 
about star- formation rates and metal enrichment, without 
a star f ormation hiatus. We also r epeat the warning pro- 
vided in lSchonrich fc Binnevl l|2009l ): it is important to bear 
in mind that our simulation data is a kinematically unbi- 
ased sample, while most similar observational plots are for 
samples that are kinemat i cally b iased in favour of thick-disc 
stars. In fact. iBovv et al.l l|2012l ) shows that the bi-modality 
in [a/Fe] is greatly diminished once the populations are cor- 
rectly weighted. 

We next analysed the properties of the two populations 
of stars as defined by the aforementioned linear separation. 
In Fig [3] we plot the star formation history of stars in the 
solar annulus (black line); the formation history of stars 
above the dashed line in Figure [5] (the putative thick disc) 
are plotted in red while the stars below the dashed line in 
Figure [2] (thin disc) are plotted in blue. A simple cut at 
the low metallicity end of the thick disc trajectory through 
this abundance space also allows us to classify stars with 
[Fe/H] < — 1.3 as belonging to the stellar halo, although we 
will see that this population may be better described as in- 
cluding a metal weak thick disc (MWTD), particularly in 
the range — 1.5< [Fe/H]< — 1.3. The star formation history 
of halo stars is plotted in yellow in Figure [3] The separation 
of the ages of stars, as selected by abundances, is clear. 

The thin disc stars are exclusively young, mostly less 
than 6.6 Gyr old, but with a small number of stars as old 
as 8.5 Gyr. Thick disc stars are mostly older than 6.6 Gyr, 
and range from 5.5 to 10 Gyr. The ages of our selected halo 
population ranges from 9.5 to 12 Gyr. The mass of stars in 
the thin, thick and halo/MWTD components in this region 
are 22.6, 8.9, and 1.3 xlO 7 M0, respectively. Dotted lines 
indicate times chosen for a pure age classification of the 
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Figure 5. The age-velocity dispersion relation for solar annulus 
stars, showing 0~ u,v,W as green diamonds, blue crosses, and red 
asterisks, respectively. 
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Figure 6. 0~w versus radius of young stars (age<200 Myrs) in 
the disc plane (|z| <0.5kpc) at three different epochs: z=1.5 (red 
line); 2=1 (yellow line); 2=0.5 (blue line). 



components, which will be explored alongside the metallicity 
cuts in the next sub-section. 

3.2 Component Abundances 

In Figure [2] we show the the distributions of [Fe/H] and 
[O/Fe] in the solar annulus stars for the thin disc (blue line), 
thick disc (red line), and halo (yellow) populations, again us- 
ing abundance cuts (solid lines) and age cuts (dashed lines) . 
This merely confirms that the selection of stars based on 
metallicity corresponds closely to the ages of the popula- 
tions. 

3.3 Component Kinematics 

In Fig [3] we plot the distribution of the rotational velocity 
(V) of the solar annulus stars for the thin disc (blue line) 
thick disc (red line), and halo (yellow) populations, using 
abundance cuts (solid lines) and age cuts (dashed lines) as 
defined in Section [3.11 The thick disc lags the thin disc by 
~25 km/s. The V distribution of halo stars shows there is 
rotation in this low-metallicity stellar population, although 
there is a hint of more than one population. The dispersion 
in U (radial velocity), V, and W (the velocity perpendicular 
to the disc plane) of the stars in the solar annulus for the thin 
disc, thick disc and halo, are (in km/s) (7tj,v,w= (41,31,25), 
(49,44,35) and (62,46,45), respectively. 

3.4 Evolution of Velocity Dispersion 

The age-velocity dispersion relation for solar annulus stars 
is plotted in Fig [S] showing <7u,V,W as green diamonds, 
blue crosses, and red asterisks, respectively. A gradual heat- 
ing is apparent to a lookba ck time of ~11 Gyr, m ore in 
line with the conc l usions of iNordstrom et al.1 j2004h than 
iQuillen fc Garnettl (|200ll ). although this is a complex issue 
that we have exp lored earlier in the context of more massive 
disc simulations (|House et al.ll201ll ). In this latter work, we 
showed that numerical heating is not causing the gradual 
heating that we are seeing in these simulations. 

We examine the evolution of velocity dispersion in more 
detail in Fig [6] where we plot 0~ w versus radius of young 
stars (age< 200 Myrs) at three different times, firstly at 



z = 1.5 (red line), a time corresponding to the early period 
of thick disc star formation, then at z = 1 (yellow line) which 
is at the end of the thick disc star formation, and finally a 
later time (z=0.5 blue line) during the thin disc formation 
epoch. At each epoch, we use young stars in the disc plane 
(|z| < 0.5 kpc) of the galaxy. It is evident that the disc stars 
that were forming at high redshift were born kinematically 
hotter than than those born at later times. Combined with 
FigQ] we can say that the stars forming at the end of the gas 
rich merger epoch create a relatively short and kinematically 
hot disc structure. We note that the extent of this old disc is 
also reflected in the z—1.5 dispersion-radius plots, with the 
jump in dispersion at R^ y ~7.5kpc, while the disc extends 
beyond 9 kpc by z=l (and beyond 12 kpc at z—Q). Yet the 
difference in dispersion at time of birth is not enough to ex- 
plain the dispersion-age relation in Figure [5] The dispersion 
radius relation is important in this respect. We will see that 
old stars in the solar annulus that have large dispersions in 
Figure [5] were preferentially born at lower radii, and hence 
with relatively high dispersion. 



4 ORIGIN OF THE COMPONENTS 

33.1xl0 7 M o of stars are in the solar annulus as defined as 
the region 7 < R xy < 8 kpc and \z\ < 0.5 kpc. Based on our 
abundance cuts, 22.6x10 7 Mq (68%) are classified as thin 
disc, 8.9xl0 7 M Q (27%) as thick disc and 1.3xlO 7 M (4%) 
as halo/MWTD. We now trace the gas from which stars in 
the solar annulus were formed. We identify five categories 
(Table [TJ for the origin of the stars, based on where they 
were formed and the source of the gas from which they were 
formed. Firstly, stars formed in proto-galaxies other than 
the central galaxy and then accreted are classified as stel- 
lar accretion (ST ACC). These accreted stars comprise just 
0.5% of stars in the solar annulus, and all were classified as 
halo stars due to their abundances being [Fe/H]< — 1.3. So 
17% of the stellar halo in the solar annulus was accreted 
in this manner. This would appear to be a fairly low frac- 
tion, and we have already stated that the halo as defined 
may include a metal weak thick disc. We make a stricter 
criteria and check retrograde halo stars, and find that still 
only 21% are direct stellar accretion. We note that accreted 
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MASS 
(10 7 M Q ) 


ST 

ACC 


CLPY 
GAS 
ACC 


I-SU 
ST 
MGE 


I-SU 

GAS 
MGE 


SMTH 
ACC 


Total 


33.1 


0.7% 


18% 


4.% 


33% 


44% 


Thin 


22.6 


0% 


16% 


0% 


28% 


56% 


Thick 


8.9 


0% 


23% 


9% 


49% 


19% 


Halo 


1.3 


17% 


28% 


47% 


7% 


1% 



Redshift 
0.5 1 



0.20 c: 



Table 1. The origin of stars in the solar neighbourhood. Stel- 
lar accretion (ST ACC): stars accreted from merged non-central 
progenitors/satellites. Clumpy Gas Accretion (CLPY GAS ACC): 
gas accreted from merged non-central progenitors/satellites. In- 
situ Star Merger (I-SU ST MGE): stars born in central galaxy 
during the merger epoch. In-situ Gas Merger (I-SU GAS MGE): 
gas in central galaxy during the merger epoch. Smooth Accretion 
(SMTH ACC): gas accreted subsequent to the merger epoch. 



stars have metallicities [Fe/H]< 1.8, so we finally make an 
even more strict criteria by limiting our sample to those with 
[Fe/H]< 1.8 and find that 31% of stars with these low metal- 
licities are directly accreted. No disc stars (thick or thin) in 
our simulation have been accreted directly. 

Stars forming from gas which is in proto-galaxies other 
than the central galaxy when it is accreted, but that does 
not form stars until after the accretion, are categorised as 
clumpy gas accretion (CLPY GAS ACC). 28% of halo stars, 
17% of thick disc stars, 10% of thin disc stars, and 13% of all 
stars in the solar annulus were formed via this mechanism. 

Stars which formed from gas that was originally ac- 
creted smoothly but form stars in the central galaxy during 
the "merger epoch" , were classified as in-situ stellar merger 
(I-SU ST MGE) stars: 47% of halo stars and 9% of thick disc 
stars are classified as due to an in situ star merger. Stars 
which form from gas that is originally accreted as smoothly 
to the central galaxy, but do not form stars until subsequent 
to the "merger epoch", are classified as in situ gas merger 
(I-SU GAS MGE) stars: 7% of halo stars, 49% of thick disc 
stars, and 28% of thin disc stars are in situ gas merger stars, 
with such stars comprising a total of 33% of all solar annulus 
stars. The remaining stars form from gas that is smoothly 
accreted to the central galaxy and are classified as being due 
to smooth accretion (SMTH ACC). 

The halo component is dominated by a combination of 
direct stellar accretion, and stars which form in the central 
galaxy (in situ) prior to the end of the "merger epoch" and 
are "knocked" into the halo by the merger events. A signif- 
icant number of the latter halo stars formed from gas that 
was originally accreted as gas from non-central progenitors. 
The remaining 7% of the halo stars, as selected by their low 
metallicity, formed from gas which was in the central galaxy 
during the gas-rich mergers but formed stars during the tur- 
bulent period which followed the merger epoch. These stars 
actually have a net average rotation of 30 km/s, and may 
be associated better with a metal weak thick disc. 

What we call 'thick disc stars' are those which formed 
predominantly from gas involved in gas-rich mergers, with 
49% being in the central galaxy prior to the end of the 
merger epoch and 17% in the smaller proto-galaxies ac- 
creted during this merger epoch. However, the thick disc 
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Figure 7. The star formation history of all stars of the simulated 
galaxy is shown as the dashed black line, plotting the fraction of 
stars born per Gyr in bins of 0.5 Gyr. Overplotted is the star 
formation history of stars in the disc (black line). We then show 
the star formation history of disc stars according to different for- 
mation radii: 3<R,F orm <5 (red), 5<Rp orm <7 (light blue), and 
7<R.Form<9 (blue). 



stars which form from this gas do so subsequent to this 
merger epoch, while the galaxy is 'settling'. Only ~20% 
of the thick disc stars in the simulation are formed from 
smoothly accreted, non-merger, gas - i.e., gas accreted to 
the galaxy that is not involved in merger events. No thick 
disc stars are directly accreted during mergers. Our findings 
highlight that gas-rich mergers are a dominant process in 
the formation of thick disc populations. Note that this is 
not at odds with studies which claim that cold flow accre- 
tion is the dominant source of a galaxies baryons. Rather, 
these studies, as well as ours, need to be interpreted care- 
fully: baryons accr eted prior to a major merger in studi es of 
cold accretion (e.g. lKeres et al.l2005l : lBrooks et"al 1l2009l ) arc 
counted as "smooth accretion" , which is dominated by cold 
flows at high redshift, and at all times that the virial mass 
is less than 10 M©. Yet, significant amounts of this gas 
which is classified as cold flow accretion can be involved in 
significant, and even major, gas-rich mergers subsequent to 
being accreted. The importance of gas-rich mergers is thus 
severely underestimated if these studies are superficially in- 
terpreted. Of course this goes both ways: the 49% of the 
thick disc that we have classified as "in situ gas mergers" 
were formed form gas that first accreted to the disc region 
smoothly. 

Thin disc stars are dominated by smoothly accreted gas. 
We note that a not insignificant amount of gas that feeds the 
thin disc does come from gas-rich mergers. Much of this is 
recycled to the disc via the hot halo, after being ejected from 
the star forming regions of the galaxy during starbursts. This 
large scale galactic fountain process allows the recycled gas 
to gain angular momentum (B12) and aids in the suppres- 
sion of the ubiquitous G-dwarf problem (Pilkington et al. 
2012 in prep). 
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Total 


a-Old 


Metal Poor 
a- Young 


Metal Rich 
a- Young 


hi 


2.73 


2.31 


4.07 


2.74 


h z 


0.7 


0.9 


0.7 


0.5 



Table 2. The scale-length (h;) and scale-height (h z ) of stars as 
selected by abundances, defining a-Old as those above the dashed 
line in Fig. 3, and a- Young as those below the line, a- Young are 
then split using a simple metallicity cut at [Fe/H]=— 0.68. Metal 
Poor a- Young have [Fc/H]< —0.68 and Metal Rich as a- Young 
have [Fe/H]> —0.68. Stars are from the region 4 < H X y < 12kpc 
and < |z| < 4kpc, while scale heights use stars in the annulus 
7 < R X y < 8kpc and are measured out to 4kpc. 

4.1 Thick Disc to Thin Disc Transition 

There is no significant "event" occurring ~6— 7 Gyr ago that 
we can associate with the separation of the thin and thick 
discs. However, the galaxy's star formation rate was rela- 
tively high between ~7— 10 Gyr ago, and decreased some- 
what near the time of the transition to thin disc star for- 
mation. In Figure [7] we plot the star formation history of all 
stars of the simulated galaxy as the black dashed line, plot- 
ting the fraction of stars born per Gyr in bins of 0.5 Gyr. We 
also show the star formation history of stars in the disc re- 
gion (7 < R xy < 8 kpc and \z\ < 0.5 kpc, black line), and the 
star formation history of disc stars according to different for- 
mation radii : 3<RForm<5 kpc (red), 5<RForm<7kpc (light 
blue), and 7<R,F rm<9kpc (blue). RForm, calculated using 
physical units, is the distance from the centre of the galaxy 
at which the star forms. This demonstrates the inside-out 
nature of the disc growth, which will have consequences for 
the abundances that we discuss below. The red line of Fig- 
ure[7]is also interesting in this respect. The stars in the solar 
annulus at z=0 that formed in the region 3<RForm<5kpc 
(the region where thick disc stars predominantly formed) 
also has a star formation rate that declined around the time 
of transition to thin disc star formation. Stars formed in 
regions further out have more a gradual build-up of star 
formation, meaning a longer enrichment timescale. 



5 STRUCTURE OF ABUNDANCE SELECTED 
POPULATIONS 

Here we mimic I Bow et al.l (|201 lh . and analyse the scale- 
length (h;) and scale-height (h z ) of stars as selected by abun- 
dances. We use our thick and thin disc separation as shown 
in our Fig. 3 to delineate the a-old population as those 
above the dashed line, and a-young as th ose below (using the 
nomenculture used in I Bow et al. l JpOllh . a-young are then 
split using a simple metallicity cu 10 at [Fe/H]=-0.68. Metal 
poor a-young (MPa-Y) have [Fe/H]< —0.68 and metal rich 
a-young (MRa-Y) have [Fe/H]> —0.68. In this section we 
have extended the region of stars beyond the solar annu- 
lus in order to measure these properties, but the metallicity 
cuts remain valid separators of the populations. Stars taken 
are from the region 4 < R xy < 12 kpc and < |zj < 4 kpc, 

2 results are not sensitive to this exact value 
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Figure 8. The evolution through abundance space of stars in 
different radial regions. The symbols represent median values of 
different regions in the disk, with R X y (kpc) 4-6 kpc ('+' sign), 
6-8 kpc (star), 8-10 kpc (diamond), 10-12 kpc (triangle). At any 
given time, the inner region is more chemically evolved than the 
outer region. This explains the longer scale length of the less 
evolved metal poor a-young population (MPa-Y) compared with 
the metal rich a-young population (MRa-Y). 



while scale heights use stars in the annulus 7 < R^ < 8 kpc 
and are measured out to 4 kpc. The scale-lengths and scale- 
heights of the three populations are listed in Table 2. 

As already noted, our a-old (thick disc) population has 
a relatively short scale-length (2.31 kpc) and large scale- 
height (0.9 kpc) compared to the total population which has 
values of 2.73 kpc and 0.7 kpc respectively. The features of 
the a-young are a little l ess simple, and int erestingly match 
very well the findings in iBow et all i|201ll ). Firstly looking 
looking at scale-heights within the solar annulus, the MPa-Y 
has a significantly larger scale-height than the MRa-Y pop- 
ulation. This simply reflects their average age, with metal 
rich stars being generally younger. Yet the MPa-Y popula- 
tion has a significantly larger scale-length than the MRa-Y 
population. This is the result of the difference in chemical 
abundance timescales between the inner and outer disks, re- 
lating to inside out disc growth. The inner regions are the 
most chemically evolved, meaning that there remains greater 
relative numbers of un-evolved MPa-Y stars out further in 
the disk which results in their population having a longer 
scale length. Figure [8] highlights this. The symbols represent 
median values of different regions in the disk, with R xy 4- 
6kpc (+ sign), 6-8 kpc (star), 8-10kpc (diamond), 10-12 kpc 
(triangle). At any given time, the inner region is more chem- 
ically evolved than regions further out. 

If migration is efficient, stars migrating to the outer disk 
will pollute the gas there with SN la, i.e. r educe the O/Fe 
ratio in the outer regions (see Figure 14 in lLoebman et all 
l201ll ). Thus, the degree of migration may be refleceted in 
the observed difference in scale-length between MPa-Y and 
MRa-Y populations. 
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Figure 9. The evolution of the distribution of formation radii for 
stars that end up in the solar annulus: 7 < R xy < 8, \z\ < 0.5. 

6 MIGRATION 

Migration proves to be another key process affecting the 
distribution of galactic components. Figure [5] shows where 
stars formed that end up in the stellar annulus. The stel- 
lar annulus is demarked by the horizontal dashed lines. The 
dotted vertical lines show the divisions between age defined 
thin disc, thick disc and halo populations. The mean forma- 
tion radius for thick disk stars is closer to the centre of the 
galaxy than the thin disc stars. Clearly, thick disc stars are 
born closer to the centre of the galaxy than the later form- 
ing thin disc stars. We do not further examine the mode of 
migration here - i.e., we do not distinguish the contributions 
of churning, blurring, scattering from the bar, or co-rotation 
scattering in thi s study. Briefly, we can compare our plot and 
Figures 3 fc 5 of lLo"ebman et al.l (|201lh : the degree of migra- 
tion is larger than in their Figure 5, the case where little 
co-rotation scattering occurs due to weak spiral structure, 
but is lower than in Figure 3, where significant co-rotation 
scattering occurs. 

Figure [TD] shows a component by component analysis of 
migration. For this analysis, we revert to selecting the com- 
ponents based on metallicity. The solid black line shows the 
total distribution of formation radius. The red line shows 
the thick disk and the blue line shows the thin disk distri- 
butions. 

We separate outlier stars out of the thick and thin disk 
populations, to try to understand why some young stars 
share abundances with the thick disc population, while some 
old stars share abundances with the thin disc. The young 
(age < 5.2 Gyr) thick disk population is shown as the dot- 
ted red line, and the old (age > 7.8 Gyr) thin disc population 
is shown as the dotted blue line. Clearly, old thin disc stars 
were born closer to the centre of the galaxy than the gen- 
eral disc population, and at even smaller radii than the thick 
disc population. This is explained by the fact that the inner 
region is more che mically evolved becau se the disc grows 
inside-out (see also IPilkington et al.ll2012T ). so the [a/Fe] ra- 
tio becomes lower at an earlier time in the innermost region 
than in regions further out. Thus, old stars born in this in- 
ner region have thin disc abundances. Contrastingly, young 
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Figure 10. The thin black line shows the distribution of radii of 
star formation of solar annulus stars, with the distribution of thin 
(blue dotted) and thick (red dotted) overplotted. The blue and 
red lines show the extremely old (age>7.8Gyr) thin and young 
(age <5.2 Gyr) thick disc stars. 

thick disc stars (based on abundances) migrated from outer 
regions, which were relatively un- evol ved chemically. This 
is sim ilar to findings in the models of dSchonrich fc Binnevl 
l2009h 

We also examine the dispersion of stars in these pop- 
ulations, and find that, for the old thin disc <7\y=31 km/s, 
comparable to other stars of similar age. aw of young thick 
disc stars (migrated from outer regions) is ~36 km/s, sim- 
ilar to the rest of the thick disc population. We compare 
with young stars that have migrated from the inner regions 
(RForm<5 kpc - i.e., the region from which thick disc stars 
predominantly migrated) and find <tw=25 km/s. This indi- 
cates that migration alone will not build a thick disc in our 
simulation, and highlights the importance that the original 
old disc is born relatively hot, as shown in Figure [6] 



7 DISCUSSION 

7.1 Milky Way Comparisons 

We do not want to over-emphasise any comparison to the 
Milky Way that can be drawn from our study, and re- 
mind the reader again that the simulated galaxy is sig- 
nificantly lower mass, and will have a different merger 
history. We do want to emphasise that this simulated 
galaxy has a low mass, low metallicity, stellar halo. 
The majority of previous simulations that have analysed 
the origins of stell ar components have had overly mas- 
sive spheroids (e.g. lAbadi et al.l 120031; IZolotov et ai.l [20091 : 
iKobavashi fc Nakasatdl201ll : iTissera et al.ll201lr ) (from the 
accretion of an excess of stars) and suffered from an inabil- 
ity t o match the empirical stell ar mass-halo mass relation 
(e.g. IPiontek fc Steinmetzl koiln . The mass and metallic- 
ity of the stellar halo p rovides strict const r aints on galaxy 
form ation simulations l|Brook et alj |2004a: Oka moto et alj 
2005). We showed that we were able separate thin and thick 
disc st ars based on abundan ces (and age), in a similar man- 
ner to iNavarro et~ail <|201lf > for the Milky Way. The qual- 
itative similarities between the simulation and the Milky 
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Way of having two tracks for the two populations through 
the [O/Fe], [Fe/H] plane are indeed interesting, and we at- 
tribute the change from thick disc to thin disc star formation 
to a drop in the star formation ra te, associated with the end 
of the rapid gas accretion phase (|Keres et al.ll2005h and the 
end of the gas rich merger epoch. The degree of separation of 
these tracks may be an instructive diagnostic: the two tracks 
in the simulation are closer together (separated by ~0.01- 
0.05 dex in [O/Fe]) t han the Milky Way tracks (separated 
by ~0.2 dex in [O /Fe] iBensbv et atbood ; iReddv et al1l2006l : 
lFuhrmannll2008l ). which has been interpreted by chemical 
evolution models as indicating that t he Milky Way experi- 
enced a hiatus in star f o rmation (e.g. Chiappini et al.ll200ll ; 
iFenner fe Gibsonlliool ). ISchonrich fe Binnevl l |2009l ) attain 
bimodal distributions in [O/Fe] as a consequence of their as- 
sumptions about star- formation rates and metal enrichment, 
without a star formation hiatus (their star formation is as- 
sumed to exponentially decay). Yet we remind the reader 
that the degree of offset in the observed relations in Milky 
Way may, at least partially, be due to the kinematic selection 
which is biased toward thick-disc stars. 



7.2 Components 

We have made a component analysis of a simulated disc 
galaxy. We categorise stars as belonging to the thin disc, 
thick disc, and halo by simple cuts in abundances of [O/Fe] 
and [Fe/ H] in a manne r simila r to that done for observed 
stars in iNavarro et all (|201lh . In qualitative sense, the 
properties of the components are reasonable analogues of 
the Milky Way components, in terms of velocity dispersion, 
rotation velocities, and metallicy. We trace the origin of gas 
from which stars in the different components form: 

Halo: In the solar neighbourhood, most stars form in 
situ but are knocked into the halo during the merger 
epoch. This is quantitatively high er than the fraction of in 
situ halo stars as suggested in (jZolotov et al.l l 2009). who 
predicted a significant fraction of halo stars are formed in 
situ, but found that accretion remains the dominant source 
of halo stars. We attribute this difference to the low stellar 
mass fractions of the accreted protogalaxies/satellites in 
our simulations, which better match empirical re lations 
between stellar and halo mass (|Brook et al.ll2012bl ). This 
also results in the low mass, low-metallicity halo of the 
current simulation. By c ontrast, using physi cal recipes 
similar to those used in IZolotov et al.l (|2009l ) has been 
shown to result in significantly too many stars compared to 
the empirical relatio n between stellar and total halo masses 
|Sawala et al.|[201ll ) in low mass galaxies, exaggerating the 
number of accreted stars and resulting in spheroids that 
are more massive and more metal-rich than observed in 
the Milky Way. In our simulated galaxy there is still a 
significant fraction (~17-31%) of local halo stars that are 
directly accreted from smaller proto-galaxies/satellites. 

Thick Disc: Thick disc stars all form in situ. Again this is in 
contrast to numerous studies of simulated Milky Way mass 
galaxies, which claim significant numbers of accreted thick 
disc s tars (e.g. Abadi et al.l 120031 ; iKobavashi fe Nakasatd 
l201ll : iTissera et al.l 201 if ). We note that those previous 
studies have stellar spheroids that are significantly more 



massive than observed in the Milky Way. Significant 
numbers of accreted thick disc stars have not been shown 
to form withou t correspond i ng ma s sive spheroids also b eing 
accreted (e.g. lAbadi et all 120061 : ITissera et"aH l201ll ). In 
our simulated galaxy, around 77% of thick disc stars form 
from gas that is directly accreted to the central galaxy, 
while 23% is accreted as clumpy gas - i.e., gas which is 
part of an accreted proto-galaxy/satellite. Yet, most of 
the thick disc forming gas is accreted prior to the end of 
the 'merger epoch', meaning that most (~80%) of thick 
disc gas is involved in a gas-rich merger at high redshift. 
Thick disc stars form as the gas settles into a disc at the 
end of the merger epoch. The thick disc has a relatively 
short radial scalelength and is kinematically 'hot' as it forms. 

Thin Disc: The thin disc stars form primarily from gas 
that is smoothly accreted to the central galaxy subsequent 
to the end of the merger epoch. A significant number of 
thin disc stars form from gas accreted during the merger 
epoch, either directly to the central galaxy or as part 
of an accreted proto-galaxy, which is cycled through the 
galaxy's hot coron a in a large scale galactic fountain (see 
iBrook et aT]|2012al '). 

Transitions: In our simulation, the end of the merger 
epoch marks the transition from halo star formation to 
thick disc star formation. No "event" marks the transition 
from thick to thin disc in our simulation. We do note that 
the thick disc formation time is longer in our simulation 
than most estimates o f the Milky Wa y thick disc (^3.5 Gyr 



co mpared to 1-2 G vrs iFuhrmann e.g 



yr_ 

e.g. I2009I . although IBensbv et al 



20081 : ISmilianic et all 



2007 claims a timescale of 



3 Gyr). This may be due to the lower mass of our simulated 
galaxy, which would result in a longer time for the thick 
disc to settle down. This would be consistent with obser- 
vations which indicate that t hick discs are more prominen t 
in lower mass disc galaxies (jYoachim fc Dalcantorl l2006h . 
Alternatively, it may of course be a problem with our model. 

Inside-out Chemical Evolution: The thick disc (with 
an a-old population, i.e. stars above the dashed line in Fig. 
3) has a scale length of 2.31 kpc at z = 0. The inside-out 
growth of the thin disc (a-young stars, below the dashed 
line in Fig. 3) results in differences in chemical abundance 
timescales between the inner and outer disks. The inner 
regions are the most chemically evolved, meaning that 
there remains greater relative numbers of un-evolved metal 
poor a-young stars out further in the disk, which results 
in that population having a longer scale length (4.07 kpc) 
compared to the metal-rich a-young stars (2.74 kpc). 

Migration: The stars classified as thick disc in our model 
preferentially were born closer to the centre of the galaxy, 
and migrated to the solar neighbourhood. Our model is fully 
cosmological, yet it is chosen to have a quiescent merger pe- 
riod at late (z < 1.5) times. At t hese t imes, it thus resembles 
the models of lLoebman et all (|201lT ). and it is interesting 
that their collapse model also initially forms a relatively 
short, kinematically hot, disc at high redshift before the 
thin disc formation begins. The dissipative nature of the gas 
rich merger epoch in our case and the collapse of a rotating 
sphere in their case (and in the study of Sam land fc Gerhardl 
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l2003h has resulted in information of the early differences be- 
tween the mod e ls bei ng lost, at least in kinematics alone. 
lLoebman et al.l l|201ll ) also found that simple cuts in abun- 
dance space could result in populations with thick and thin 
disc properties. Detailed chemical "fingerprints" that will 
become available in surveys such as GALAH may com- 
bine with GAIA kinematic data to unravel star forma- 
tion c ites (jFreeman fc Bland-Hawthornll2002l ; lDe Silva et all 
2006). We also note th at the initial conditions imp osed in 
the dynamical model of Sc honrich fc Binnevl l|2009l ) consist 
of a relatively hot and short disc. Our simulations suggest 
that migration is part of the story of thick disc formation, 
but not necessa rily to the exclusion of o ther theories. The 
recent paper of l|Liu fc van de Venll2012h finds evidence for 
both migration and gas rich mergers in the relationships be- 
tween eccentricities and abundances. 
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